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Abstract. The purpose of this paper is to present simple and fast methods for computing 
control points for polynomial curves and polynomial surfaces given explicitly in terms of 
polynomials (written as sums of monomials). We give recurrence formulae w.r.t. arbitrary 
affine frames. As a corollary, it is amusing that we can also give closed-form expressions in 
the case of the frame (r, s) for curves, and the frame ((1, 0, 0), (0, 1, 0), (0, 0, 1) for surfaces. 
Our methods have the same low polynomial (time and space) complexity as the other best 
known algorithms, and are very easy to implement. 



1 



1 Introduction 



Polynomial curves and surfaces are used extensively in geometric modeling and computer 
aided geometric design (CAGD) in particular (see Ramshaw fOj, Farin [SIE], Hoscliek and 
Lasser 7], or Piegl and Tiller 9j). One of the main reasons why polynomial curves and 
surfaces are used so extensively in CAGD, is that there is a very powerful and versatile 
algorithm to recursively approximate a curve or a surface using repeated affine interpolation, 
the de Casteljau algorithm. However, the de Casteljau algorithm applies to curves and 
surfaces only if they are defined in terms of control points. There are situations where a 
curve or a surface is defined explicitly in terms of polynomials. For example, the following 
polynomials define a surface known as the Enneper surface: 



Thus, the problem of computing control points from polynomials (defined as sums of mono- 
mials) arises. If control points can be computed quickly from polynomials, all the tools 
available in CAGD for drawing curves and surfaces can be applied. This could be very use- 
ful in problems where a curve of a surface is obtained analytically in terms of polynomials 
or rational functions, for example, problems involving generalizations of Voronoi diagrams, 
or motion planning problems. When dealing which such problems, it is often necessary to 
decide whether curves segments or surface patches intersect or not. As is well known (for 
example, see jSl IS), there are effective methods based on subdivision (exploiting the fact 
that a Bezier curve or surface patch is contained within the convex hull of its control points) 
for deciding whether Bezier curve segments or surface patches intersect. Simple and fast 
methods for computing control points might also be also useful to teach say. Math students, 
to learn computational tools for drawing interesting curves and surfaces. Now, it turns out 
that the problem of computing control points can be viewed as a change of polynomial ba- 
sis, more specifically as a change of basis from the monomial basis to bases of Bernstein 
polynomials. Algorithms for performing such changes of basis have been given by Piegl 
and Tiller |9j. More general algorithms for performing changes of bases between progressive 
bases and Polya bases are presented in Goldman and Barry j3] and Lodha and Goldman 
(H]. These algorithms compute certain triangles or tetrahedras whose nodes are labeled with 
certain multisets, and are generalizations of the de Casteljau and the de Boor algorithm. 
In this paper, we present alternate and more direct methods for computing control points 
from polynomial definitions (in monomial form) that run in the same low time complexity 
as the above algorithms (O(m^) for curves of degree m, 0{p^q^) for rectangular surfaces of 
bidegree {p,q), and 0{m^) for triangular surfaces of total degree m). Our algorithms are 
not as general as those of Goldman and Barry [S] and Lodha and Goldman jH] , but they are 
more direct and very easy to implement. 
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The paper is organized as follows. In section |2l we review briefly the relationship between 
polynomial definitions and control points. We begin with the polarization of polynomials in 
one or two variables, and then we show how polynomial curves and surfaces are completely 
determined by sets of control points. In the case of surfaces, depending on the mode of 
polarization, we get two kinds of surfaces, bipolynomial surfaces (or rectangular patches) 
and total degree surfaces (or triangular patches). Efficient methods for computing control 
points are given in the next three section: polynomial curves in section |3J bipolynomial 
surfaces in section HI and polynomial total degree surfaces in section El Some examples are 
given in section IHl 



2 Control Points 
2.1 Polynomial Curves 

The deep reason why polynomial curves and surfaces can be handled in terms of control 
points is that polynomials in one or several variables can be polarized. This means that 
every polynomial function arises from a unique symmetric multiaffine map. A detailed 
treatment of this approach can be found in Ramshaw ^Oj, Farin 0l2J, Hoschek and Lasser 
[7j, or Gallier We simply review what is needed to explain our algorithms. 

Recall that a map /: M'' — > is ajfine if 

/((l_A)a + A6) = (l-A)/(a) + A/(6), 
for all a,be M'^, and all A G M. A map /: M^>< - ■ ■><R^ ^ is multiaffine if it is affine in 

m 

each of its arguments, and a map f:W^x ■ ■ ■ x R'^ ^ is symmetric if it does not depend 

m 

on the order of its arguments, i.e., /(a7r(i), . • . , an{m)) = f{(^i, ■ ■ ■ , Q-m) for all ai, . . . , a^, and 
all permutations vr. We also say that a map f:MF x M'^ — > M'^ is {p,q) -symmetric if it is 
symmetric separately in its first p arguments and in its last q arguments. 

Let us first treat the case of polynomials in one variable, which corresponds to the case 
of curves. Given a (plane) polynomial curve F: M ^ of degree m, 

x(t) =Fi(t), 

yit) = F2(t), 

where Fi{t) and F2{t) are polynomials of degree < m, it turns out that F: R — >■ comes 
from a unique symmetric multiaffine map /: — R^, the polar form of F, such that 

F{t) = f{t,. . . ,t), foralHGR. 
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Furthermore, given any interval (r, s) (affine frame), the map /: M"^ — >■ is determined 
by the sequence (60, • • • , &m) of m + 1 control points 

m—i i 

where < i < m. Using hnearity, in order to polarize a polynomial of one variable t, it is 

^ j terms in the sum 

E (n'<)- 

/C{l,...,m} iei 

\I\=k 

the polar form f^{ti, . . . , tm) of the monomial t'^ with respect to the degree m (where k <m) 
is given by 

I 1 /C{l,...,m} ie/ 

V A; y l^l=fc 

2.2 Polynomial Surfaces Polarization 

Given a polynomial surface F: — > R^, there are two natural ways to polarize the polyno- 
mials defining F. 

The first way is to polarize separately in u and v. If p is the highest degree in u and q is 
the highest degree in we get a unique (p, g)-symmetric degree (p + q) multiaffine map 

such that 

F{u,v) = f{u,...,u;v,...,v). 



We get what is traditionally called a tensor product surface, or as we prefer to call it, a 
bipolynomial surface of bidegree {p, q) (or a rectangular surface patch). 

The second way to polarize is to treat the variables u and v as a whole. This way, if F 
is a polynomial surface such that the maximum total degree of the monomials is m, we get 
a unique symmetric degree m multiaffine map 

/: (R^)"* ^ R^ 

such that 

F(u,v) ^ f((u,v),...,(u,v)). 



We get what is called a total degree surface (or a triangular surface patch). 

Using linearity, it is clear that all we have to do is to polarize a monomial w^'f 

It is easily verified that the unique {p, g)-symmetric multiaffine polar form of degree p + g' 

fP'l{ui,...,Up;vi,...,Vg) 

of the monomial u'^v'' is given by 



^1(^1 /c{i,...,p},|/|=h \ia / 
h \ k jc{i,...,,},\j\=k 



The denominator ( f I ( ? 1 is the number of terms in the above sum. 



h J \ k 

It is also easily verified that the unique symmetric multiaffine polar form of degree m 

f]^kiiui,Vl),...,iUm,Vm)) 

of the monomial u'^v'^ is given by 



m \ f m — h . ^, , 

h I \ k I \i\=h.lJ\=k,inJ-- - 



The denominator ^'^^^"^^^^"^^^^^^^^ A;))^^ number of terms in 
the above sum. 

2.3 Control Points For Polynomial Surfaces 

Let = {{i, j, k) e \ i + j + k = m}. Given an affine frame Arst in the plane (where 
r,s,te are affinely independent points), a polynomial surface F: — > R^ of total degree 
m specified by the symmetric multiaffine map 

/: (M^)™ ^ m3 

is completely determined by the family of points in R^ 

i j k 

where {i,j,k) e A^. 
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These points are called control points, and the family {bij^k \ G Am} is called a 

triangular control net. 

Let (ri,si) and (^2,52) be any two affine frames for the affine line M. A bipolynomial 
surface F: — of bidegree (p, q) specified by the (p, g)-symmetric multiaffine map 



called a rectangular control net. 

Thus, to compute control points, in principle, we need to compute the polar forms of 
polynomials. However, this method requires polarization, which is very expensive. In the 
following sections, we give recurrence formulae for computing control points efficiently. As 
a corollary, in the case of any affine frame (r, s) or of the affine frame ((1,0), (0, 1), (0,0)), 
it is possible to give closed-form formulae for calculating control points in terms of binomial 
coefficients. 

3 Computing Control Points For Curves 

We saw in section |21 that the polar form of a monomial t^ with respect to the degree m is 



is completely determined by the family of {p + l)(g + 1) points in 




where < i < p and < j < g. 

These points are called control points, and the family {b. 



< i < p, < j < q} is 




Letting a] 
equations: 



.m 



( 



m 



k 



) 



/™, it is easily verified that we have the following recurrence 




.m— 1 



m— 1 



if 1 < < m; 

if = and m > 0; 

otherwise. 



The above formulae can be used to compute inductively the polar values 
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tions, we observed that computing and dividing by I , ) is faster than computing /™ 



The computation is reminiscent of the Pascal triangle. Alternatively, we can compute 
/™ directly using the recurrence formula 

i-m _ ~ k) i-m-l 

Ik ~ Jk ' ''"i Jk-1 ' 

m m 

where 1 < k < m. When writing computer programs implementing these recurrence equa- 

m 
k 

directly using the above formula. This is because the second method requires more divisions. 

Given {ti, . . . ,tm), computing all the scaled polar values cr^(ti, . . . ,U), where I < k < 
i and 1 < i < m, requires time 0{m?). The naive method using polarization requires 
computing Y^^q 2* = 2™"*"^ — 1 terms. To compute the coordinates of control points, we 
simply combine the a^. Specifically, the coordinate value of the control point hj contributed 
by the polynomial XlILo ^^^^ (where m < m) is 

-1 



k—O m—j j 



where (r, s) is an affine frame. Given (ti, . . . , tm), our algorithm computes the table of values 
al(ti, . . . ,ti), where 1 < k < i and 1 < i < m, and thus, it is very cheap to compute these 
sums. 

Remark: Given a polynomial curve F of degree m specified by the sequence of control 
points {Bq, . . . ,bm) over (0,1), it is well known (see Farin |Sl Hoschek and Lasser [7j, 
or Piegl and Tiller [HI) that F{t) can be expressed in terms of the Bernstein polynomials 

ITT' \ ■\in-k+k 



F{t) = B^{t)h, + --- + BZ{t)h^. 

It is also well known that the Bernstein polynomials B^{t), . . ., B^{t) form a basis of the 
vector space of polynomials of degree < m. Thus, it is also possible to compute the control 
points of F by expressing the polynomials involved in the explicit polynomial definition of 
F in term of the basis Bernstein polynomials. Such algorithms were given by Goldman and 
Barry [Sj. Our algorithm has the same complexity and is more direct. 

It is also easy to derive closed- form formulae for any affine frame (r, s). 

Theorem 3.1 // the total degree is m and there are r occurrences where t = p and s occur- 
rences where t = q, then 



fk 



k 



Proof . It can be shown by induction using the above recurrence equations. □ 
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4 Computing Rectangular Control Nets 

As we saw in section |2l the polar form of the monomial u^v'^ with respect to (p, q) is 



fh'iiui, ...,up;vi,..., = - — ^ — ^ ^ f n f n 



p \ q 



^3 



/c{i,...,p},|/|=h Vie/ / Vjgj 
h I \ k I jc{i,...,5},|j|=fc 



Letting o"^'^ ~ /i ^ ^ easily verified that we have the following recurrence 



equations: 



a 



h,k 



(^h,k'' ^ + Up^h-\l ^ + ^9<,fc-'i ^ + «P^</^L};Li if 1 < < P and 1 < < g 
(Tq + '"g'^o fc-i if /i = < p and 1 < /c < g 



+ '^p^h-io li 1 <h <p and /c = < g 

1 ' ' if /i = = 0, p > 0, and g > 

^ otherwise. 



Observe that the recurrence formula is a sort of generalization of the Pascal triangle. 
Alternatively, prove that /^'| can be computed directly using the recurrence formula 

fP,q _ jP ~ h){q — k) p_i^g_i h{q — k) j,p-l,q-l , jP ~ rp-l,q-l I ^ rp-l,q-l 

Jh,k- Jh,k + '^pJh-i,k + '"qJh,k-i + p^'^p'"<iJh-i,k-v 

where 1 < h < p and 1 < /c < g, 



where = < p and 1 < < g, and 

/h,0 - p Jhfl + p"p//i-l,0' 

where 1 < h < p and A; = < g. As in section |3J we found that computing cr^'^ and dividing 
by ^ ^ ^ ^ ^ ^ is faster than computing /^'^ directly. 

Given {ui, . . . ,Up;vi, . . . ,Vg), using the recurrence equations, computing all the scaled 
polar values 

where l<h<i,l<k<j,l<i<p, and 1 < j < g, can be done in time 0{p^q'^). The 
naive method using polarization requires computing XliLo =o 2*^"' ~ 2^"'''''*'^ — 1 terms. 



To compute the coordinates of control points, we combine the cr^'l- Specifically, the 
coordinate value of the control point bij contributed by the polynomial X^ILo YlT=o ^h,kU^v^ 
(where m < p and m < g) is 




h=0 k=0 



where (ri,si) and (^2,52) are affine frames. Our algorithm computes the table of values 
(J™^™'(mi, . . . , Urn] ^i, . . . , Vm) , and thus, it is very cheap to compute these sums. 

When the affine frames (0, 1) are used, the following theorem gives closed-form formulae 
for the polar values with respect to the bidegree (p, q). 

Theorem 4.1 // there are s occurrences where u = 0, r occurrences where v = 0, and all 
the other occurrences of u and v have the value 1, then 



p — s \ 1 q — r 
h [ k 



p \ I q 
h \ k 



Proof . It can be shown by induction using the above recurrence equations. □ 

It is well known that F can be expressed in terms of control points and (products of) 
Bernstein polynomials (see Farin 0121; Hoschek and Lasser j7j, or Piegl and Tiller As 
in the case of curves, the methods of Lodha and Goldman [H] have the same complexity as 
ours, but our method is more direct. 

5 Computing Triangular Control Nets 

As we saw in section |2l the polar form of the monomial u^v^ with respect to the total degree 



m IS 



fj!^kiiuuvi),---,ium,vj) = - — ^ — V fn^' I fn 



^3 



/ujc{i,...,m} Vie/ / VjeJ 

h I \ k I \i\=h,\j\=k,inj=$ 



^ j ( ^ j fj^f,, it is easily verified that we have the following 



recurrence equations: 



+ w-^rTfc + ^mc^rfc-i if Z^, A; > O and l < + A; < m. 



o'h.k = ii h = k = and m > 

otherwise. 
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The above formulae can be used to compute inductively the polar values 



m \ I m — h 
h [ k 



The computation consists in building a tetrahedron of values reminiscent of the Pascal tri- 
angle (but 3-dimensional) . Alternatively, we can compute fj^f, directly using the recurrence 
formula 

rm _ — h — k) nm-l I _^ nm-l 

Jh,k- ^ Jh,k ^'^rnJh-l,k^ ^'"mJh,k-V 

where h,k > and 1 < h + k < m. Again, we found that computing cr^j. and dividing by 
'^^^"^^^^is faster than computing /^j. directly. 

Given {{ui,vi), . . . , {um,Vm)), using the recurrence equations, computing all the scaled 
polar values 

where h,k > 0, 1 < h + k < i, and 1 < i < m, can be done in time O(m^). The naive 
method using polarization requires computing J2^o 3* = (3"*"*'''^ — 1)/2 terms. To compute the 
coordinates of control points, we simply combine the crj^^j^.. Specifically, the coordinate value 
of the control point bij^k (where i+j+k — m) contributed by the polynomial X^/i+Km ^h,iu'^'^^ 
(where m < m) is 




where d = m — k — l and r = (ri, r2), s = (si, S2) and t = (ti, ^2) are the vertices of the affine 
frame. Our algorithm computes the table of values cr^ ^((iii, Vi), . . . , (uj, Vj)), and thus, it is 
very cheap to compute these sums. 

When the affine frame ((1,0), (0,1), (0,0)) is used, the following theorem gives closed- 
form formulae for the polar values with respect to the total degree m. 

Theorem 5.1 Assume that m — r + s + t, with r occurrences of (1,0), s occurrences of 
(0,1), andt occurrences o/(0, 0) (and no occurrences of (1,1)). Then 



fJ^kii'^i^'"l)^---A'^m,Vm)) = 



r \ / s 
h [ k 



m \ 1 m — h 
h )\ k 
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Proof. It can be shown by induction using the above recurrence equations. □ 

As in the previous case, it is well known that F can be expressed in terms of control 
points and (trivariate) Bernstein polynomials (see Farin O Hoschek and Lasser [7j, or 
Piegl and Tiller |9|). The methods of Lodha and Goldman jH] have the same complexity as 
ours, but our method is more direct and very easy to implement. 



6 Examples 

We wrote an implementation in Mathematica of a program computing polar values for curves, 
using the recurrence equations of section IHl It works for an arbitrary affine frame (r, s) . 

As an example, consider the curve of degree 10 given by 

_ 4t(l-t^)^(l-14t' + t^) 

_ 8t^(l-t=^)(3-10t^ + 3t^) 
^~ (l + t2)5 • 

Using the above program, the following control polygon w.r.t. [0, 1] is obtained: 

rcpoly = {{0, 0, 1}, {2/5, 0, 1}, {18/25, 12/25, 10/9}, {1/2, 6/5, 4/3}, 
{-14/45, 71/45, 12/7}, {-45/37, 45/37, 148/63}, {-71/45, 14/45, 24/7}, 
{-6/5, -1/2, 16/3}, {-12/25, -18/25, 80/9}, {0, -2/5, 16}, {0, 0, 32}}; 

Note that the control points also contain weights, since there are denominators (see 
Ramshaw p^]], Farin [31 Ej, Hoschek and Lasser Jl^, or Gallier |3]). Here is the rational 
curve. 
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Figure 1: A rose 



We also wrote an implementation in Mathematica of a program computing polar values 
for triangular surface patches, using the recurrence equations of section This algorithm 
works for any affine frame (r, s, t). 

In Hilbert and Cohn-Vossen [Gj (and also do Carmo PP), an interesting map 7i from 
to is defined as 

(x, y, z) ^ {xy, yz, xz, x^ - y^). 

This map has the remarkable property that when restricted to the sphere S"^, we have 
z^ = 7Y(x', y', z') iff (x', y', z') = {x, y, z) or (x', y', z') = (— x, —y, —z). In other words, 
the inverse image of every point in Ti^S"^) consists of two antipodal points. Thus, the map Ti 
induces an injective map from the projective plane onto Ti^S"^), which is obviously continuous, 
and since the projective plane is compact, it is a homeomorphism. Thus, the map Ti allows 
us to realize concretely the projective plane in M^, by choosing any parameterization of the 
sphere S*^, and applying the map Ti to it. For example, the following parametric definition 
specifies the entire projective plane over [—1, 1] x [—1, 1]: 

16'ut'^(l — u^) 
_ Suv{u^ + l){v'^ - I) 

^~ (m2 + 1)2(^2 + 1)2 ' 
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_ Av{l - u^){v^ - 1) 
_ ^v^u" - + 1) 

Using our algorithm, tiie following net of degree 8 over the affine frame ((1, 0, 0), (0, 1, 0), 
(0, 0, 1)) is obtained: 

projSnet = 

{{0, 0, 0, 0, 1}, {0, 0, -1/2, 0, 1}, {0, 0, -14/15, 2/15, 15/14}, 

{0, 0, -20/17, 6/17, 17/14}, {0, 0, -120/101, 60/101, 101/70}, 

{0, 0, -1, 4/5, 25/14}, {0, 0, -11/16, 15/16, 16/7}, {0, 0, -1/3, 1, 3}, 

{0, 0, 0, 1, 4}, {0, 0, 0, 0, 1}, {0, -1/7, -1/2, 0, 1}, 

{4/45, -4/15, -14/15, 2/15, 15/14}, {4/17, -28/85, -20/17, 6/17, 17/14}, 

{40/101, -32/101, -120/101, 60/101, 101/70}, 

{8/15, -6/25, -1, 4/5, 25/14}, {5/8, -1/8, -11/16, 15/16, 16/7}, 

{2/3, 0, -1/3, 1, 3}, {0, 0, 0, 0, 15/14}, {0, -4/15, -7/15, 0, 15/14}, 

{20/121, -60/121, -105/121, 9/121, 121/105}, 

{10/23, -14/23, -25/23, 9/46, 46/35}, 

{240/331, -192/331, -360/331, 108/331, 331/210}, 

{80/83, -36/83, -75/83, 36/83, 83/42}, {10/9, -2/9, -11/18, 1/2, 18/7}, 
{0, 0, 0, 0, 17/14}, {0, -32/85, -7/17, 0, 17/14}, 

{9/46, -16/23, -35/46, -1/46, 46/35}, 
{27/53, -89/106, -50/53, -3/53, 53/35}, 
{36/43, -100/129, -40/43, -4/43, 129/70}, 

{12/11, -6/11, -25/33, -4/33, 33/14}, {0, 0, 0, 0, 101/70}, 

{0, -48/101, -34/101, 0, 101/70}, 

{56/331, -288/331, -204/331, -40/331, 331/210}, 

{56/129, -44/43, -98/129, -40/129, 129/70}, 

{16/23, -144/161, -120/161, -80/161, 23/10}, {0, 0, 0, 0, 25/14}, 

{0, -14/25, -6/25, 0, 25/14}, {8/83, -84/83, -36/83, -16/83, 83/42}, 

{8/33, -38/33, -6/11, -16/33, 33/14}, {0, 0, 0, 0, 16/7}, 

{0, -5/8, -1/8, 0, 16/7}, {0, -10/9, -2/9, -2/9, 18/7}, {0, 0, 0, 0, 3}, 

{0, -2/3, 0, 0, 3}, {0, 0, 0, 0, 4}}; 

Note that the control points also contain weights, since there are denominators. If we 
project the real projective plane onto a hyperplane in R^, either from a center or parallel to 
a direction, we can see a "3D shadow" of the real projective plane in R^. For example, one 
of the projections is the cross-cap, whose control net is 

projnet4 = 

{{0, 0, 0, 1}, {0, 0, 0, 1}, {0, 0, 2/15, 15/14}, {0, 0, 6/17, 17/14}, 
{0, 0, 60/101, 101/70}, {0, 0, 4/5, 25/14}, {0, 0, 15/16, 16/7}, 
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{0, 0, 1, 3}, {0, 0, 1, 4}, {0, 0, 0, 1}, {0, -1/7, 0, 1}, 

{4/45, -4/15, 2/15, 15/14}, {4/17, -28/85, 6/17, 17/14}, 

{40/101, -32/101, 60/101, 101/70}, {8/15, -6/25, 4/5, 25/14}, 

{5/8, -1/8, 15/16, 16/7}, {2/3, 0, 1, 3}, {0, 0, 0, 15/14}, 

{0, -4/15, 0, 15/14}, {20/121, -60/121, 9/121, 121/105}, 

{10/23, -14/23, 9/46, 46/35}, {240/331, -192/331, 108/331, 331/210}, 

{80/83, -36/83, 36/83, 83/42}, {10/9, -2/9, 1/2, 18/7}, {0, 0, 0, 17/14}, 

{0, -32/85, 0, 17/14}, {9/46, -16/23, -1/46, 46/35}, 

{27/53, -89/106, -3/53, 53/35}, {36/43, -100/129, -4/43, 129/70}, 

{12/11, -6/11, -4/33, 33/14}, {0, 0, 0, 101/70}, {0, -48/101, 0, 101/70}, 

{56/331, -288/331, -40/331, 331/210}, {56/129, -44/43, -40/129, 129/70}, 

{16/23, -144/161, -80/161, 23/10}, {0, 0, 0, 25/14}, 

{0, -14/25, 0, 25/14}, {8/83, -84/83, -16/83, 83/42}, 

{8/33, -38/33, -16/33, 33/14}, {0, 0, 0, 16/7}, {0, -5/8, 0, 16/7}, 

{0, -10/9, -2/9, 18/7}, {0, 0, 0, 3}, {0, -2/3, 0, 3}, {0, 0, 0, 4}} 

The cross-cap is shown below. 



1 




Figure 2: The cross-cap surface 



Acknowledgement. I wish to thank Deepak Tolani for very helpful comments, in 



14 



particular on the use of Bernstein polynomials. 
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